Study of the topological Hall effect on simple models 
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Recently, a chirality-driven contribution to the anomalous Hall effect has been found that is 
induced by the Berry phase and does not directly involve spin-orbit coupling. In this paper, we will 
investigate this effect numerically in a two-dimensional electron gas with a simple magnetic texture 
model. Both the adiabatic and non-adiabatic regimes are studied, including the effect of disorder. By 
studying the transition between both regimes the discussion about the correct adiabaticity criterium 
in the diffusive limit is clarified. 



I. INTRODUCTION 

When studying the Hall effect in fcrromagnets, an 
anomalous contribution was found giving rise to a non- 
zero Hall effect even in the absence of an externally ap- 
plied magnetic field. Spin-orbit coupling was invoked 
to explain this effect, as it gives rise to two scattering 
mechanisms (skew scattering^'^ and side jump'^) leading 
to a preferential scattering direction that is different for 
spin-up and spin-down. Since the spin subbands in fcr- 
romagnets are unequally populated, this spin-dependent 
scattering can indeed give rise to a charge Hall effect 
in fcrromagnets, while in normal semiconductors these 
mechanisms explain the much-discussed spin Hall effect"*. 

Quite recently, one has found that an anomalous Hall 
effect can exist even in the absence of spin-orbit inter- 
action in some classes of frustrated ferromagnets like 
pyrochlore-type compounds with non-coplanar magnetic 
moments in the elementary cell^. In these systems, the 
Hall effect is a result of the Berry phased an electron 
acquires when moving in a system with nontrivial (chi- 
ral) spin texture. Since the effect can be attributed to 
the topology of the magnetization field, the term topo- 
logical Hall effect was coined^. Theoretical calculations 
on this effecliSiS have mainly concentrated on the adi- 
abatic regime, where the electron spin aligns perfectly 
with the local magnetization during its movement. In 
this regime, by doing a transformation that aligns the 
quantization axis with the local magnetization direction 
at every point in space, the original Hamiltonian can be 
mapped onto a model of spinless electrons moving in an 
effective vector potential2i^. However, only a few papers 
have dealt with the non-adiabatic limit (using perturba- 
tion theory^") and even less is known about the transition 
between the two regimes. 

In this paper, we will show numerical calculations on 
the topological Hall effect in a two-dimensional electron 
gas (2DEG). It was proposed that in this system an arti- 
ficial chirality can be created e.g. by the stray field of a 
lattice of ferromagnetic nanocylinders placed above the 
2DEG, having the advantage that all relevant parameters 
are controllable up to some extenliX. In the current paper 
however, we will only consider some very simple model 
chiralities in the system. Nevertheless, with these models 
we are able to investigate numerically the existence of the 



topological Hall effect. More importantly, the transition 
point between the non-adiabatic and adiabatic regime 
will be studied for different values of the elastic mean 
free path. In the diffusive regime, it is found that the 
transition occurs slower as we increase the amount of 
disorder. This sheds some light on a long-standing dis- 
cussion about the relevant adiabaticity criterium; it is in 
agreement with the one derived by Stern* ^ and speaks 
against the criterium put forward by Lossi^. 

The paper is subdivided as follows. In the next section, 
we will first give a more detailed description of the topo- 
logical Hall effect, discussing the transformation of the 
original Hamiltonian in the adiabatic limit onto a model 
of spinless electrons moving in an effective vector poten- 
tial. After different adiabaticity criteria are presented 
in section UTTl a discussion concerning the calculation of 
Hall resistances and resistivities follows. In particular, a 
special phase averaging method is introduced for describ- 
ing large (i.e. system length L ^ phase coherence length 
L^) disordered systems; technical details of this method 
are given in the appendix. The main results of the paper 
can be found in section Ivl 



II. TOPOLOGICAL HALL EFFECT 

In order to explain what the topological Hall effect is 
about, we will consider a two-dimensional electron gas 
(2DEG), subjected to a spatially varying magnetization 
M(r). The Hamiltonian of such a system can be written 



H = - — V^-ga.Mirl 



(1) 



with m* the effective electron mass, a the vector of Pauli 
spin matrices, and g a coupling constant. We assume that 
the amplitude of the magnetization is constant; only its 
direction is position dependent: M(r) — Mn(r), with n 
a vector of unit amplitude. 

For our numerical purposes, it is necessary to obtain 
a tight-binding description of this system by discretizing 
the Schrodinger equation on a square lattice with lattice 
constant a. The tight-binding equivalent of Eq. Q reads 
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where i,j label the lattice sites, a,/? are spin indices, 

is the hopping amplitude, and n^ — n(rj). The 



t 



2ma'^ 



first summation runs over nearest neighbors. 

When the coupling g is large enough, the adiabatic 
limit is reached meaning that spin flip transitions will be 
absent and the electron spin will stay aligned to the local 
magnetization direction. In that case, the solutions for 
the spin-up spinors (with respect to the local magnetiza- 
tion) are given by 



coi,(4.)exp(-ife) 
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with {9i,4>i) the spherical coordinates of the local mag- 
netization direction n^. By projecting the Hamiltonian 
onto the subspace spanned by these spinors, one can de- 
rive an effective HamiltonianSiii 
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where the hopping amplitudes are now 
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In this expresssion, Qj is the angle between the vectors 
TO-i and Wj. The angle 7^ needs some more explanation. 
Suppose an electron makes a closed trajectory around a 
lattice cell like depicted in Fig. ^ Since its spin follows 
adiabatically the local magnetization direction, the elec- 
tron will pick up a Berry phase that is equal to half the 
solid angle f2 subtended by the magnetization directions 
at the four corners of the cell&, and the quantity ^ij is 
related to this Berry phase. 

To make things more clear, let's consider a flux 



$/$o = ^ 



(6) 



piercing through a lattice cell and forget about the mag- 
netization for a moment ($0 = h/e is the magnetic 
flux quantum). The influence of such a fluxtube can 
be described in a tight-binding model by a Peierls sub- 
stitutionji^ which would change the hopping phase as 
t — > texp(— ie/?i/ A • dl), with the integral evaluated 
along the hopping path. This is similar to the phase fac- 
tor in Eq. ((SJ with 
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meaning that in principle the Berry phase effect can be 
described by an effective vector potential. E.g., one could 
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FIG. 1: Mapping of the tight-binding Hamiltonian for an 
electron in a magnetic texture to a spinless electron moving 
around a fluxtube. Grey color of the vertices on the right 
picture stand for the change in hopping amplitude, while the 
arrows denote a change in hopping phase due to the fluxtube. 



make a choice of gauge for A such that the hopping 
phase on all the vertices above the fluxtube change as 
t -^ iexp(i27r$/$o) = texp(if7/2) (see Fig. [2). This 
gauge choice would correspond to the Landau gauge in 
a situation where the effective fluxes through all lattice 
cells would be equal (i.e. there is a homogeneous field). 
In short, the effect of the 7^ can thus be described by 
having a distribution of effective magnetic fluxes piercing 
through the lattice cells. As such, for a given gauge, the 
quantities 7ij are defined in a unique way. 

In summary, the effect of the transformation on the 
tight-binding Hamiltonian will be twofold: first, the hop- 
ping amplitude between neighbors changes, depending on 
the angle between the magnetization directions, and sec- 
ondly, there is a Berry phase effect that effectively can 
be described in terms of a fluxtube distribution, where 
the value of the flux (in units $0) through a single lattice 
cell is given by ^, and 11 is the solid angle subtended by 
the magnetization directions at the four corners of the 
lattice cell. 

By means of this mapping it is clear now that the mag- 
netization texture can indeed give rise to a Hall effect. 
It should be noted that the effective flux is given by a 
solid angle, so that the topological Hall effect can only be 
nonzero for non-coplanar textures. Furthermore, there is 
no need to invoke any kind of spin-orbit coupling for the 
topological Hall effect to appear. 



III. ADIABATICITY CRITERIUM 

In order to be able to do the transformation between 
the Hamiltonian describing an electron moving in a mag- 
netization texture [Eq. (^, which will be referred to as 
the magnetization model from now on], and the effec- 
tive Hamiltonian in Eq. Q) [called the flux model from 
now on], the electron spin has to stay aligned to the lo- 
cal magnetization direction during its movement. This 
adiabatic regime will be reached when the spin preces- 
sion frequency lOs is large compared to the inverse of a 
timescale r that is characteristic for the rate at which the 



electron sees the magnetization direction change. 

If transport through the nanostructure is balhstic (i.e. 
there is no disorder), then this timescale is given by r = 
(./vp, where vp is the Fermi velocity of the electrons, and 
^ is the length scale over which the magnetization changes 
its direction substantially. In this case, the adiabaticity 
criterium is thus 



Vp 



>1. 



(8) 



When introducing disorder in the system, it is clear that 
this criterium is still valid as long as the mean free path 
Im is larger than ^. However, when going to the strongly 
diffusive regime /,„ < ^, two different timescales r ap- 
pear in the literature and there is still a discussion going 
on about the relevant one.^iii24iSii& In analytical calcula- 
tions on rings, Stern has found that the relevant timescale 
would be the elastic scattering time t = Tm, leading to 
what is referred to as the pessimistic criterium in the 
literatureii 



Q»r-. 

t'Tn. 



(9) 



However, other papers by Loss and coworkersiSi have put 
forward that the relevant timescale is the Thouless time 
'''th = {^l^mY '''rn, i-S. the time the particle needs to dif- 
fusive through a distance ^. This would lead to a less 
stringent adiabaticity criterium since rth > Tm in the dif- 
fusive regime, and it is therefore known as the optimistic 
criterium: 



Q»y. 



(10) 



In a paper of van Langen et al. the pessimistic criterium 
was confirmed by a semi-classical analysis^^, and later 
numerically by Popp et al}^ In section ^ we will have 
a closer look at the transition point between the adia- 
batic and non-adiabatic regimes ourselves and find also 
confirmation for the criterium in Eq. ^. 



IV. CALCULATION OF THE HALL 
RESISTANCE 

In order to observe the topological Hall effect, we will 
numerically calculate the Hall resistance in a geometry 
like shown in Fig. [3 The Hall resistance Rh in this 
geometry is defined as 



1 

Rh — ^(^12,34 — ^34,12)j 

where we use the common notation 

Vk-Vi 
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for a measurement where current is supplied through con- 
tacts i and j, and the voltage difference V^ — Vi is mea- 
sured, fixing Ik = Ii = 0. Making the difference between 



3® 



©2 



II ® I • I * t * II 

3^©— »®^©— »(• 

] ® I « I « I ® I 

•)—►©<—©—»©<— 6 

1 « t • i • t * 1 

0<^©^©<^©^0 . 

1 ® || • I ® t ® I ® II ® 4 

13— >©<—©->©<— 6i 

FIG. 2: Geometry of a single section used to calculate the 
Hall resistance. A possible magnetization texture is shown. 
In the leads, the magnetization is chosen to point out of the 
plane of the paper. 



these two four-point resistances is equivalent to defining 
the Hall resistivity pn = \{Pxy — Pyx) as the part of the 
resistivity tensor that is anti-symmetric with respect to 
time reversalii. 

Within the Landauer-Biittiker (LB) formalism, the 
four-point resistances in Eq. (|12|l can be expressed as& 



R 



■ij^kl 



h TkiTij — TiiTf 
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(13) 



as a function of transmission probabilities Tpq between 
the leads. In this expression, D is an arbitrary 3x3 
subdeterminant of the matrix A relating the currents and 
the voltages on the leads: Ip — ^ ApqVq. 

In principle, one could proceed now by calculating Rh 
as a function of the adiabaticity parameter Q and com- 
pare the rate of convergence for samples with different 
mean free paths in order to find the relevant adiabaticity 
criterium in the diffusive regime. However, since we make 
use of the LB-formalism, it is assumed that the whole 
structure is phase coherent. It is well known that the 
resistances in such structures are quantitatively strongly 
dependent on the actual configuration of the impurities in 
the system even when they are characterized by the same 
mean free path. Although this is an integral part of the 
physics of mesoscopic systems, it makes a quantitative 
comparison of Rh between samples with different mean 
free paths useless. We would like to compare properties 
of a macroscopic system, i.e. a system with a finite phase 
coherence length in which such fluctuations are absent. 

Therefore, some kind of (phase) averaging over differ- 
ent disorder configurations has to be introduced to find 
a description of the transport properties in terms of a 
macroscopic material constant, like the Hall resistivity. 
Some care should be taken in defining such an averaging 
procedure; e.g., just calculating the mathematical aver- 
age Rh = jj J2i=i ^H ^^ ^^^ Hall resistances i?|f found 
for N different impurity configurations does not give a 
quantity that is directly related to the Hall resistivity 
PH of a macroscopic system. 

The basic idea behind the averaging procedure we have 
chosen is that a macroscopic system (L ^ L^) can be 
thought to consist of smaller phase coherent sections of 
size L « L^. For every smaller section, we can use the 
LB-formula to derive its transport properties, and the 
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FIG. 3: Basic cell for the texture of the magnetization (a). 
Every arrow corresponds to the magnetization direction at a 
single lattice site. By scaling up once, a more smoothly vary- 
ing magnetization texture is obtained, resulting in smaller 
solid angles subtended by a single lattice cell (b). In (a) all 
arrows lie in the plane of the paper, while in (b) the mag- 
netization is shown projected onto this plane: smaller arrows 
mean a shorter projection. 
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FIG. 4: Hall resistance as a function of energy for a single 
section with magnetization texture obtained by scaling up the 
cell in Fig. EK twice. 



properties of a macroscopic system can then be found 
by attaching such sections in an incoherent way. In our 
case, the smaller sections from which a macroscopic sys- 
tem will be built up are structures with the geometry 
in Fig. [2 Since a more detailed discussion of our par- 
ticular averaging procedure and the corresponding Hall 
resistivity pu is rather technical, it is given in a separate 
appendix at the end of the paper. 



V. RESULTS 

Before discussing the main results of the paper, a few 
words should be said about the model we have chosen for 
the magnetic texture in our numerical calculations. We 
started from a 4 x 4 cell like depicted in Fig. |2^. Only 
the direction of the magnetization changes from site to 
site like depicted in the figure, while its magnitude stays 
constant. However, when mapping this texture to a flux 
model (like described in section ^J ^ one finds that the 
4x4 cell comprises a total flux of 4 $07 corresponding to 
an effective flux per lattice cell of $ = l/4$0j which is 
quite large. One could decrease this value by decreasing 
the solid angle subtended by the magnetization direction 
vectors at the corners of each lattice cell, since the flux 
per lattice cell is proportional to it [see Eq. ((HJ] . In order 
to do so one can scale up the cell of Fig.|3^ by introducing 
extra lattice sites in between the original sites, and by 
interpolating the magnetization direction at the new sites 
between the directions of the neighboring sites. Scaling 
up the cell of Fig.|2Ji once, one then finds a magnetization 
cell of 8 X 8 sites, like depicted in Fig. |3j3. Of course the 
total flux comprised by this cell is still 4 $0; but this time 
distributed over 64 sites, which corresponds to an average 
flux per lattice cell of 1/16 <i>o. This magnetizion cell 
can then in turn again be scaled up; every time we scale 
up the cell, it will comprise 4 times the original number 



of sites so that the average flux per lattice cell will be 
decreased by four. It should be noted that although the 
effective flux per lattice cell for the 4x4 magnetization 
unit cell we started from is homogeneously distributed 
(exactly 1/4 $0 per lattice cell), this is not anymore the 
case for the cells found by scaling up the flrst one. 

As a flrst step, we will have a look at a small sections 
with the geometry in Fig. [3 In the rest of the paper, 
we will model these sections by a tight-binding lattice 
comprising 65 x 65 sites. The leads connected to them 
have a width of 30 sites. The magnetization in the leads 
connected to the sample is chosen to point out of the 
plane of the paper on every site. When translating this 
to the flux model, it means that there is zero magnetic 
flux per lattice cell in the leads. In the device itself, a 
magnetization cell of 16 x 16 sites is used, found from 
scaling up the original cell in Fig. 0^ two times. We 
can fit 16 such magnetization cells in one of our 65 x 65 
samples. 

In Fig. 21 we have calculated the Hall resistance Rh as 
defined in Eq. Hill) for a single structure like depicted in 
Fig. [21 and for energies of the incoming electrons ranging 
from up to 2 i above the bottom of the spin up band. 
The spin splitting due to the magnetization was chosen as 
large as A = gM ~ 100 1, making sure that we are in the 
adiabatic regime (no disorder is present in the system). 
The results clearly show a nonzero Hall resistance Rh\ 
in fact, one can clearly observe the integer quantum Hall 
effect. This indeed proves that a Hall effect is present 
although we do not have any magnetic flux present, nor 
any kind of spin-orbit coupling. The topological Hall 
effect seen here is purely due to the Berry phase an elec- 
tron picks up when moving through a structure with a 
magnetization texture with nonzero spin chirality. 

As explained in section^ one can map the model with 
the magnetization texture to a spinless flux model. For 
the magnetization texture we have chosen, the effective 
average flux per lattice cell will be <I>/$o — 1/64. This 
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FIG. 5; Hall resistivity as a function of the adiabaticity factor Q for different values of the mean free path. 



corresponds to a cyclotron radius of 10 a (for the Fermi 
energy Ep = li), which is much smaller than the sam- 
ple size so that the integer quantum Hall effect could 
indeed be expected. When doing the mapping onto the 
flux model by calculating the effective hopping parame- 
ters like in Eq. Q for our particular choice of magneti- 
zation texture, and then calculating the Hall resistance 
Ru , we found indeed perfect overlapping between the re- 
sistances of the two models. It can be claimed that the 
perfect overlap between the two models does not come 
as a surprise because in the quantum Hall regime the 
Hall resistance is quantized into exact plateaus. How- 
ever, calculations were done even outside the quantum 
Hall regime, and the overlap between the two models is 
always exact within our numerical accuracy. Further- 
more, the main point here is the fact that a Hall effect 
can be observed that is due to the electron adiabatically 
following a certain magnetization structure, and this is 
clearly demonstrated by Fig. ^ 

Now that the topological Hall effect has been demon- 
strated numerically, we can proceed to study the transi- 
tion from the non-adiabatic to the adiabatic regime. This 
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will be done as a function of disorder strength (mean free 
path) , in order to check which of the two adiabaticity cri- 
teria in section IhH are correct. As disorder is introduced 
in the system (within the Anderson model) , a phase aver- 
aging procedure like set out in the previous section (and 
described fully in the appendix) has to be done in or- 
der to get quantitative information. For every separate 
value of the mean free path, we have calculated the trans- 
mission coefficients for 500 small sections (Fig. ^ with 
different impurity configurations. Then we have wired 
together 4900 sections randomly chosen from these 500 
in a square 70 x 70 array, and calculated the resistivity 
of the obtained structure like described in the appendix 
(Fig. [Tji). For the magnetization texture, we have cho- 
sen a 64 X 64 cell found by scaling up the cell in Fig. [Jt- 
three times. A single section thus comprises a single cell 
of the magnetization. The shortest distance over which 
the magnetization changes its direction by an angle tt is 
given by ^ « 22 a for this particular cell. The Fermi en- 
ergy was fixed to Ep = It above the bottom of the spin 
up subband. 
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Fig. Eh- shows the Hall resistivity pH defined in 
Eq. I|A4I) for mean free paths ^m = 32 a and Im — 64 a, as 
a function of the adiabaticity parameter Q = ojs^/vp = 
hA£^/vp. For both values of the mean free path the 
adiabatic limit is reached simultaneously for values of 
Q w 20, where the Hall resistivity becomes independent 
of Q. This is in good agreement with the adiabaticity 
criterium in Eq.|Hl Furthermore, it was checked that the 
adiabatic value of the Hall resistivity coincides with the 



one calculated by mapping the magnetization model to 
the flux model. It should also be noted that pn is inde- 
pendent of the mean free path {pn « 5 x 10^'^ -^) which 
is essentially what one might expect from a simple Drude 
model. Another feature on the figure is that pn does not 
increase linearly for small Q; instead it stays very small 
up to some point where it changes abruptly. This point 
coincides with a value for the spin splitting A = Ep. 
Thus adiabaticity is reached on a rather short scale as 
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FIG. 6: Hall resistivity as a function of the adiabaticity factor 
Q in the intermediate regime. Decreasing the mean free path 
leads to a slower convergence to the adiabatic value. All plots 
are shifted so as to converge to 0. 



soon as the Fermi energy lies below the spin down sub- 
band. 

Figure 03 shows the same plot, but now for mean free 
paths in the diffusive regime [1^ < S,), namely /„ = 9a 
and Im = 16 a. Again, the Hall resistivity increases 
abruptly around Q « 20 to its adiabatic value. The 
resistivity pn obtains again the same value as in Fig. EJi 
around pH ^ 5 x 10"'^ h/e"^ and is thus clearly inde- 
pendent of the mean free path. Although pn changes 
abruptly, it can be seen that the real adiabatic value is 
reached slower for the lower of the two mean free paths 
(Jm = 9 a) for which pn first overshoots its adiabatic 
value, and then slowly converges to it. This difference is 
made more visible in Fig. El where we plotted the differ- 
ence between the Hall resistivity pn and the adiabatic 
value it reaches (so that all curves converge to 0), for 
mean free paths of ^m = 9 a, 12 a and 16 a. Again, it is 
clear that increasing the amount of disorder leads to a 
slower convergence to the adiabatic regime. As such, we 
confirm that the pessimistic adiabaticity criterium, which 
states that adiabaticity is reached more slowly upon de- 
creasing the ratio lm/£, [see Eq. 0], is the correct one. 
Although for our limited range of allowed parameters (we 
should have Im > Xp so that localization effects do not 
start playing a role) both adiabaticity criteria do not dif- 
fer very much quantitatively, Q S> 0.45 versus Q ^ 2 
for Im = 9 a, the optimistic criterium would predict the 
opposite behavior. 



VI. CONCLUSIONS 

In this paper, we have shown numerical calculations 
confirming the existence of a fully topological Hall ef- 
fect, that is due to the Berry phase an electron picks 
up when moving adiabatically in a non-coplanar mag- 
netization texture. In the adiabatic regime, the govern- 



ing Hamiltonian can be mapped onto a model of spinless 
electrons moving in a magnetic flux; both models indeed 
give the same results for the Hall resistance/resistivity. A 
closer look at the transition point between non-adiabatic 
and adiabatic regime revealed a rather abrupt transition 
upon increasing the exchange splitting. The transition 
takes place around the point where the spin-down sub- 
band becomes depopulated. Furthermore, we were able 
to find confirmation for the pessimistic adiabaticity cri- 
teriumii by looking at the transition point for different 
mean free paths in the strongly diffusive regime. In this 
regime, a special method for phase averaging should be 
introduced for getting rid of conductance fluctuations, 
which enables us to describe the transport properties of 
a large system (L > L^) in terms of a Hall resistivity. 



APPENDIX A: PHASE AVERAGING 
PROCEDURE 

The main idea behind our averaging procedure is that 
a macroscopic system can be built up from smaller phase 
coherent subsections that are attached in an incoherent 
way. The small sections we would like to start from were 
depicted in Fig. |21 However, there is a subtle point that 
should be considered first: it is known that at every inter- 
face between a lead and the electron reservoir feeding this 
lead, a so-called contact resistance is present (see e.giSi 
and references cited there). Essentially such contact re- 
sistance results because on one side (the reservoir) cur- 
rent is carried by an infinite number of modes, while on 
the other side (the mesoscopic lead) there are only a finite 
number of modes transporting current. When building 
up a system of macroscopic size from the smaller sections, 
such contact resistance effects in the current- voltage re- 
lations of a single section are unwanted since they are 
characteristic of the leads, and not the sample itself. 

In order to get rid of contact resistance effects, we will 
make an eight-terminal structure by attaching four extra 
voltage probes as in Fig. [7^. Since these extra voltage 
probes do not draw any current, there will be no voltage 
drop over their contact resistances, and they will measure 
the voltage at the point where they are attached. Placing 
them like in Fig. [7^, one is able to measure the voltage 
drops over the sample, excluding the voltage drops over 
the contact resistances of the original four current carry- 
ing leads (Fig. [71d). Relating now the currents through 
the original leads to the voltages measured with the volt- 
age probes, one gets rid of the contact resistance effects. 

Formally, this can be done with the Landauer-Biittiker 
formalism as follows. First we write down a set of linear 
equations relating the currents and voltages at all eight 
leads: 
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where I = (/i,...,/4) is a vector containing the cur- 
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FIG. 7: Scheme of the phase averaging procedure. Four ex- 
tra voltage probes are attached to the original four-terminal 
structure (a). At every reservoir/lead interface a contact re- 
sistance is present; these are depicted in grey in (b). Since 
the voltage probes do not draw current, they will measure 
voltage drops over the sample, not including the contact re- 
sistances. Every section can now be considered as an effective 
four-terminal box characterized by a set of equations between 
currents through the original four leads, and the voltages mea- 
sured by the extra voltage probes (c). Such boxes are wired 
up to form a macroscopic system, and its resistivities can be 
found by the van der Pauw technique (d) . 



rents through the original four terminals, and Ip = 
(/pi , . . . , Ipi) are the currents through the voltage probes; 
the same notation convention is used for the voltages on 
the leads. The 4x4 matrices A, B, C and D consist of 
transmission coefficients between all eight leads, and are 
found directly from the Landauer-Biittiker equations. 
Since the voltage probes do not draw current, we find 



Ip = = CV + i?Vp, 



(A2) 



which can be used to express the currents through the 
four current-carrying leads as a function of the voltage 
on the attached voltage probes: 



1= [B-AC-^D]Yp. 



(A3) 



Doing so, we have found current- voltage relations for the 
original four-terminal structure, but since we relate cur- 
rents through the original leads to voltages at the voltage 
probes, we have gotten rid of contact resistances. 



Subsequently, a large number of such structures (with 
different impurity configurations) are wired together as 
shown in Fig. [Tji. Every single section is treated as a 
classical four-terminal circuit obeying current-voltage re- 
lationships of the form (|A3|) [see also Fig. |7|;]. Doing so, 
we introduce an effective phase breaking event at every 
connection point between neighboring sections: only cur- 
rent and voltage information is kept at these points, and 
any phase information of electrons flowing out of the sec- 
tion is lost. In other words, an artificial phase coherence 
length L^ corresponding to the length of a single small 
section is introduced in the system. 

By applying the correct current conservation laws at 
the intersection points in Fig. [7|i, one can derive a rela- 
tionship between the currents and voltages on the four 
terminals at the corners of the connected structure (la- 
beled 1 to 4 in the Fig. [Tji). When enough sections are 
attached, this relationship will be independent of the im- 
purity configurations in the isolated sections. The prop- 
erties of the system can then be expressed in terms of re- 
sistivities by making use of the van der Pauw techniqueiS. 
First, the four-terminal resistances i?i2,34, ^34,12, ^14,23 
and -R42,3i between the corners of the large structure are 
calculated. With these, the Hall resistivity is defined 

]^ 19.20 



Ph — Rh — x(-Ri2,34 — -^34,12) 



(A4) 



while the longitudinal resistivity p^ can be found from 
solving the equationiS 

exp (-7ri?i4,23/PL) + exp (-7ri?42, 31 /pl) = 1 (A5) 

A small comment should be made here. Because of 
the attachment procedure described above, one cannot 
expect the resistivities calculated in Eqs. ljA4p and (|A5|) 
to correspond exactly to the real resistivities of a " bulk" 
system of the same size and with the same mean free 
path; in particular, the resistivities as calculated above 
depend on the width of the leads attached to the small 
sections that make up the large system, and also on the 
scheme of wiring the sections together. In theory one 
could take into account these effects (see e.g. Ref. Iia) . 
and one would find that the real "bulk" resistivity and 
the resistivity calculated above are equal up to a factor 
that is purely geometrical. Calculating this factor ex- 
plicitly is practically quite difficult; we did not proceed 
in this direction since the factor is of a purely geometric 
origin and has no physical implications. 
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